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1 Introduction 



The discovery of neutrino masses and mixings in neutrino oscillation experiments calls 
for an extension of the Standard Model (SM) of particle physics. The most elegant way 
to implement the required neutrino masses into the SM is via the seesaw mechanism pQ. 
Three right-handed neutrinos are added to the SM particle content, which, being singlets 
under the SM gauge symmetry, acquire heavy Majorana masses. By integrating out these 
right-handed Majorana neutrinos, light masses for the left-handed neutrinos are generated 
naturally. In a supersymmetric version of this mechanism, the lepton flavor violation present 
in the neutrino sector is transmitted also to the sleptons inducing charged lepton flavor 
violation (LFV). We thus focus in this work on a supersymmetric seesaw model. In the 
framework of degenerate right-handed neutrinos, LFV processes in the SUSY seesaw model 
were studied in (21 El (see also jH El El for further phenomenological studies on LFV in 
processes at low energies). 

Another particularly attractive feature of the seesaw mechanism is that it offers a natural 
possibility to generate the observed baryon to photon ratio in the universe E] , 



This process is known as leptogenesis j§] and utilizes the fact that Majorana neutrinos do not 
carry a conserved lepton number. First, a lepton number asymmetry is generated by out- 
of-equilibrium decays of the heavy right-handed Majorana neutrinos and sneutrinos (in the 
following subsumed as neutrinos) in the early universe. This lepton number asymmetry is 
later on transformed into a baryon number asymmetry via non-perturbative Standard Model 
processes induced by finite temperature transitions between topologically distinct vacua 
separated by the sphaleron barrier. If the right-handed neutrinos are generated by scattering 
in the thermal bath, leptogenesis provides a particularly straightforward prediction of the 
baryon asymmetry of the universe. The decay rates of the right-handed neutrinos and 
washout processes depend on their masses and Yukawa couplings. Hence the observed 
value of t}b (HJ may imply constraints on the seesaw model (TUl HU H2] • 
In total, the seesaw mechanism introduces 18 new parameters in addition to the Minimal 
Supersymmetric Standard Model (MSSM) parameters. A major drawback of the scheme 
is, that only 9 of these parameters describe the three masses, mixing angles and CP phases 
of the light neutrinos, which are directly accessible to the experiment. The other nine 
parameters are related to the masses and mixing of the heavy Majorana sector, which may 
be related to physics at the GUT scale. The parameters enter the renormalization group 
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evolution of the neutrino Yukawa couplings and affect the mixing of the slepton sector, 
which gives rise to observable LFV processes such as lj — > lij, with branching ratios that 
depend on the parameters of the heavy Majorana sector. 

On the other hand, these parameters enter the prediction of the baryon asymmetry via 
leptogenesis and thus also affect the lower bound on the reheating temperature of the 
universe after inflation, which governs the gravitino abundance. Since the decay of long- 
living gravitinos can dissociate light elements and thus jeopardize the successful predictions 
of Big-Bang nucleosynthesis, a low reheating temperature is desirable. 
In this work we study the constraints on the parameters of the SUSY seesaw model in a 
mSUGRA framework imposed by the combined neutrino data, successful leptogenesis with 
low reheating temperature and charged lepton flavor violation. Furthermore, we require 
that the Yukawa couplings do not become too large so that the perturbative theory holds, 
and assume hierarchical spectra for both the left- and right-handed neutrinos. 
This work is organized as follows. The theoretical framework of the SUSY seesaw model 
is outlined in Section 2. Section 3 introduces the neutrino data used in our calculations. 
Leptogenesis is explained in Section 4, followed by a discussion of the gravitino problem in 
Section 5. Constraints from leptogenesis and charged lepton flavor violation on the seesaw 
parameters are derived in Section 6 and 7. 

2 SUSY seesaw model 

The simplest supersymmetric version of the seesaw mechanism is described by the super- 
potential j3] 

W = Wmssm + \v c R T Mu c R + v c R T Y v I ■ h 2} (2) 

where Wmssm describes the MSSM, vr denotes the right-handed neutrino singlets and I the 
left-handed lepton doublets, h 2 the Higgs doublet with hypercharge +|, and Y v the matrix 
of neutrino Yukawa couplings. Below the mass scale of the lightest right-handed Majorana 
neutrino, the effective superpotential [B| is obtained by integrating out the heavy neutrino 
fields, 

W e// = Wmssm + \{Y„ I ■ h 2 f M~\Y V I ■ h 2 ). (3) 

After electroweak symmetry breaking, W e // © leads to the following mass term for the 
light neutrinos, 

m v = m^M- x m D = Yj ' M~ l Y v {hl) 2 = K(h° 2 } 2 , (4) 
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where 



v 2 sin 2 (3 is the vacuum expectation value (v.e.v.) of the relevant 



Higgs field and m^, = Y v v 2 is the neutrino Dirac mass matrix. Here it has been assumed 
that the scale of is much smaller than the scale of the heavy neutrino mass matrix, 
w,d *C M. This corresponds to the so called type I seesaw mechanism, where any direct 
(type II) contributions to the light neutrino mass matrix are neglected. Such terms can be 
generated, e.g. from v.e.v's of SU(2)-triplet Higgs fields. 

It is convenient to work in the flavor basis in which the charged lepton Yukawa matrix is 
diagonal, so that the symmetric matrix k is diagonalized by the unitary Maki-Nakagawa- 
Sakata (MNS) matrix U, 

U t kU = Diag(/«i, k 2 , k 3 ) = D K , 



the light neutrino mass eigenstates being given by 

nii = v 2 Ki, i = 1,2, 3. 
If one chooses Ki>0,U can be written as 

U = V •Diag(e-^ /2 ,e-^' /2 ,l), 



(5) 



(6) 



(7) 



where 0, 0' are Majorana phases and V can be parametrized in the standard Cabibbo- 
Kobayashi-Maskawa form, 
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cos 9 ij and Sjj = sin^-. 



From the relations Q and (0), and working in the basis where M is diagonal, M 
Diag(Mi, M 2 ,M 3 ) = D M , one obtains 



u*D K rf = {r t Jd^y v ) t {r t Jd-^y v \ 



(9) 



where R is an arbitrary complex orthogonal matrix R T R = RR T = 1. The neutrino 
Yukawa coupling matrix can thus be expressed as 



y v = Jd m rJd k iP. 



(10) 



The complex orthogonal matrix R can be parametrized by three complex angles 9^ = Xi+iyt 



R 
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(11) 



with &i = cos(^) and = sin(^). 

We focus on the simplest scenario, where the right-handed neutrinos are hierarchical 1 , 

Mi < M 2 < M 3 , (12) 

and thus the baryon asymmetry is generated essentially in the decays of the lightest right- 
handed neutrino with mass M\. Consequently, we also assume a hierarchical left-handed 
neutrino spectrum, 

mi C m 2 « m 3 , (13) 

since Yukawa couplings and hierarchical Majorana masses would have to conspire in an 
unnatural way to lead to a degenerate spectrum for the light neutrinos. 
To zero order in the small ratios M1/M3, M 2 /M 3 and m\jm 3 the Yukawa coupling (fTU|) is 
given by 

Dsinp V y 777.3 / 



with m2 — y A^?2 an d m 3 — \/^ m 23) defined through the mass squared differences mea- 
sured in neutrino oscillation experiments. Inserting the parametrization (|11|) for R one 
finally arrives at 

- — ^c 2 S i3 8u — e l<p /z V j2 + c x V* 3 . (15) 



v sin (3 



Before proceeding, we briefly summarize the renormalization group evolution of neutrino 
Yukawa couplings from m z to the GUT-scale M GUT for non-degenerate seesaw scales [T5] . 
Below Mi, the heavy Majorana neutrinos are decoupled, so that Y v = and only the effective 
light neutrino mass matrix k evolves, starting from the input value n(mz) = U*D K U\ At 
the Mi-threshold, the corresponding right-handed Majorana neutrino becomes active and 
one has 



(y„)« Mi+ = 5a UD u Ry/D K rf\ Ah _, (16) 



where Mi+ (Mi—) denotes the right (left) limit approaching the scale Mi. Furthermore, 
using (J3J, the tree- level matching condition for k at Mi is given by 

= («)« -( Y J) ttOVu ■ (17) 

v > 3 Mi+ V 111 Mx- \ v Jil Mi J Mi- V ' 

Note that the masses of heavy neutrinos also run above the respective mass thresholds. 



1 For a discussion of leptogenesis in scenarios with degenerate neutrinos see I14j . 
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Now we continue with the evolution from Mi to M 2 using the input values k\ 
and Mi|jvf 1+ . The matching 



M 2 - 



M 2 - 



M 2 - 
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(19) 



at M2 is analogous to the matching at M\. Above the M 2 threshold, the 2x2 submatrix of 
M receives small off-diagonal elements, so that at the M 3 scale M has to be diagonalized 
by an unitary matrix Um, Uj 4 M| Mg _ Um = Dm- The re-diagonalization of M leads to 
the redefinition Y v — *> U M Y U . As has been noted in j!5j . the renormalization group equa- 
tions (RGE's) are invariant under the transformations that diagonalize M. Therefore the 
matching conditions at M 3 are given by 



M 3 - 



M 3 - 



0. 



M 3 - 



ij 



M 3 - 



(20) 
(21) 



Finally, from (J2*U|) and ()2ip. we continue the evolution up to the unification scale Mx- 
With M and Y v now calculated at Mx, we now turn to the slepton sector, assuming the 
mSUGRA universality conditions 



m II, 



ITLr 



m II, 



A e = A Y, A v = A Y U , 



(22) 



at Mx, where m is the common scalar mass and A the common trilinear coupling. At 
lower scales of order of the SUSY threshold, the mass squared matrix of the charged sleptons 
has the form 

,2 



2 2 f 



ra- 



ni 



LR 
2 



(23) 



where m~ , m~ and m~ are 3x3 matrices in flavor space, m~ and m~ being hermitian. 

'l Ir Ilr 'l la 

The respective matrix elements are given by 



( m ~i R h 



[m L )ij + 5ij ( m l . + m z cos 2(3 ( -- + sin 6 W 



(m R )ij + Sij(mi. — m z cos 2/3 sin 9w), 
(A e )ijV cos/? - 5ijm h fi tan (3, 



(24) 
(25) 
(26) 



9w being the weak-mixing angle, and \i the SUSY Higgs-mixing parameter. The first terms 
on the r.h.s. of (|24|) - (|26|) receive the following contributions: 



m 2 L 



ml II + (5ml) 



+ 5m 2 L: 



(27) 
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m R = m Q l\ + (5m R ) MSSM + 6m R , (28) 
A e = AqYi + (5A e ) MSSM + 8A e . (29) 

Here, (Sm1 R ) MSSM and (SA e ) MSSM denote the usual MSSM renormalization group corrections 
[Tm iT7] which are flavor-diagonal. In addition, the right-handed neutrinos radiatively induce 
the flavor non-diagonal terms Stul^r and SA e , resulting from the RGE's 

16tt 2 ^| = mlY^Y v + Y^Y v ml + 2{Y^mlY v + m 2 h ^Y v + AlA 1 ), (30) 

16tt 2 ^ = 0, (31) 
a In /i 

l 6vr 2^£ = 2Y e Y}A u + A e YX- (32) 

For qualitative discussions, we mention here the results in the leading logarithmic approxi- 
mation, 

5ml = -^(3ml + A 2 )(YjLY u ), (33) 



■57T 2 



5m R = 0, (34) 

OA 

5A e = -j—^(YiYjLY u ), (35) 

1D7T 2 



with 

Li* = \n( 

Mi 

In the following we adopt the particle spectrum of the mSUGRA benchmark scenario SPSla 
[T8j as an illustrative example. This allows us to easily put our work in context with other 
SUSY analyses such as ^H] ■ While the numerical results depend on this choice, the generic 
features also apply to other SUSY scenarios with universality conditions such as in 



3 Neutrino masses and mixing 

An important requirement on the seesaw model is that the observed light neutrino masses 
and mixings are reproduced. This is automatically guaranteed in the R matrix parametriza- 
tion sketched in Section El where the light neutrino parameters can be used as input at the 
electroweak scale in the renormalization group evolution of the Yukawa coupling matrix. 
In the following, we use the global best fit values of the neutrino oscillation parameters 
obtained in a three neutrino framework [201 (for a more recent analysis, see j2I]) compris- 
ing KamLAND, CHOOZ, MACRO, Super-Kamiokande and SNO data as well as the first 
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spectral data from the K2K long baseline accelerator experiment. The resulting best fit val- 
ues are summarized in Tab. ^ In order to indicate the uncertainties from the experimental 
errors one will have to live with, we have used 2<r-errors expected from future measurements 
as explained in j2]. 



Parameter 


best fit 


future range 


sin 2 23 


0.52 


+0.1 

-0.1 


sin 2 6»i3 


0.005 


+0.001 

-0.005 


sin 2 0i2 


0.30 


+0.05 
-0.05 


Am 2 2 /1(T 5 eV 2 


6.9 


+0.36 
-0.36 


Am 2 3 /l(r 3 eV 2 


2.3 


+0.7 
-0.9 



Table 1: Best-fit values and 2a CL. uncertainties expected from future measurements 
of neutrino oscillation parameters. 

For the Dirac phase 8 and the two Majorana phases <fr and 0' introduced in (J7|), no experi- 
mental limits exist (compare, however [2*2"]). 

Upper bounds on the absolute mass scale of neutrinos can be obtained from tritium beta 
decay experiments, neutrinoless double beta decay searches and the neutrino hot dark 
matter effect on the cosmological large scale structure and the cosmic microwave background 
|23j . Furthermore, generation of the baryon asymmetry (Q) via thermal leptogenesis from 
the decays of the lightest Majorana neutrino Ni yields a constraint 2 on the mass scale of 
the lightest neutrino [2U : 

m l < 0.15 eV. (37) 

Moreover, as has been mentioned above, a hierarchical right-handed neutrino spectrum 
strongly favors a hierarchical left-handed neutrino spectrum. We therefore assume for con- 
sistency that next-generation double beta decay experiments will not observe neutrinoless 
double beta decay and restrict our discussion to left-handed neutrinos being lighter than 
the anticipated double beta decay sensitivity, mi < 0.03 eV. 

4 Leptogenesis 

Thermal leptogenesis as a mechanism for the generation of the baryon asymmetry of the 

universe has been discussed extensively in [26, 24J. For the sake of completeness we review in 

2 An at least mild hierarchy of right-handed neutrinos is assumed and strong phase cancellations, which 
may relax this bound, are neglected ^1] . 
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the following the most important relations and facts. The baryon to photon ratio generated 
from out-of-equilibrium decays of the heavy right-handed neutrinos is given by 



r) B ~d a Sph ei K f . (38) 

The CP asymmetry e\ generated in the decays of the lightest right-handed neutrino N\ is 
defined as |27j 

r (Ni ^h 2 + i)- rjNi -+h 2 + i) 

61 T(N 1 ^h 2 + l) + T(N 1 ^h 2 + T)' 1 j 

In the SUSY seesaw model one also has to consider the supersymmetric versions of these 
interactions involving e.g. heavy right-handed sneutrinos j22j. The CP violation in the 
decays of the Ni arises from the interference of tree-level decay diagrams with vertex and 
self-energy corrections 



* - -lj^jM^ (Y ^) f {S)' (40) 



f(x) = -fi {—^ + In i±£) . (41) 
\x — 1 X J 

M? 

For hierarchical heavy Majorana neutrinos -gjr ^> 1, j ^ 1, and apart from strong phase 
cancellations [Tij . one obtains / (^j^j — 3 If 1 , leading to ITUj 

« - J ^ (42) 

As the MNS matrix U drops out in (|42|) . it is clear that non-zero imaginary parts of the i? 
matrix elements are necessary to generate a, CP asymmetry. 
Substituting the expression (jTUj) into (|4*2|) one gets (TUj 

Using furthermore the orthogonality condition J2i^ii — 1> one approximately has [TTH 1371] 

N £ q T {m 3 - mi ). 44 

87T U| 

For a hierarchical spectrum of light neutrinos, the CP asymmetry e\ is maximal. On the 
other hand the relation (|44|) also implies a lower bound on the Mi-scale, e.g. if e\ > 10~ 6 , 
then Mi > 4 - 10 9 GeV [TO]. 
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The efficiency factor Kf takes into account the effects of the washout rate on an initial 
(B — L) asymmetry by inverse decays, and by AL = 1 and AL = 2 scattering processes. A 
reliable numerical fit for Kf for hierarchical light neutrinos in the strong washout regime is 
given by EE 121 EU E2 

/in- 2 eV\ L1±ai 
/c/(mi) ~ (1.5 ±0.7) ■ 1(T 2 f J. (45) 

The generated (B — L) asymmetry strongly depends on the effective neutrino mass 

2 • 2 a \*v*u J 11 / ac\ 

mi = v sin a . 46 

Mi 



For values of \j Am\ 2 < fh\ < J Am^ in the strong washout regime, the predictions of 
the final baryon asymmetry become independent of initial conditions and have minimal 
theoretical uncertainties [23 |2H EH] . 

Finally one also has to take into account the dilution of the asymmetry due to standard pho- 
ton production from the onset of leptogenesis until the recombination time of photons [25j, 
described by the dilution factor d = (3c/!; ec )/(4c/^ 5SM ) ~ 1/78 in the MSSM, where g^ SSM 
and g^ ec denote the degrees of freedom in the MSSM and at the time of recombination, 
respectively. 

The initial B — L asymmetry is subsequently converted to a baryon asymmetry by sphaleron 
processes. In the case of the MSSM with three fermion generations and two Higgs dou- 
blets and taking into account the chemical potentials of all particle species in the high- 
temperature phase, one obtains the conversion factor 

g 

as P h = — • (47) 

Note that this factor of roughly one third also arises in the SM with one Higgs doublet. 
Since ^ is known quite accurately and the effective neutrino mass is given by fh\ = 
J2i m i \Rii\ 5 there is a direct relation between Mi and the light neutrino masses and R 
matrix elements, 

Ml ~ Jl-M E»™*|fl"l 2 v l, (48) 

3 dasph Li mf^sm (i?fj Kf 
Note that is a direct consequence of a successful thermal leptogenesis for hierarchical 
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5 Gravitino production and reheating temperature 



A major obstacle of thermal leptogenesis can be the generation of an overabundance of 
gravitinos which causes serious cosmological difficulties [Slj. Since the couplings of the 
gravitino to ordinary matter are strongly suppressed by the gravitational scale, it has a 
very long lifetime. Nevertheless, if it is heavier than the Lightest Supersymmetric Particle 
(LSP), it can decay e.g. radiatively into a photon and a photino. These decays will occur 
after the big-bang nucleosynthesis (BBN), unless the gravitino is heavier than ~ 10 TeV j35| . 
Among the gravitino decay products are energetic photons which induce electromagnetic 
cascade processes, thereby spoiling successful BBN. Since the number of gravitinos produced 
during the reheating epoch is approximately proportional to the reheating temperature Tr, 
one can obtain upper bounds on T R depending on the gravitino mass m 3 / 2 - 
According to the analysis [HH], the upper bounds corresponding to a heavy, i.e unstable 
gravitino are given by T R < 10 7 , 10 9 and 10 12 GeV for m 3/2 = 100 GeV, 1 TeV and 3 TeV, 
respectively 3 . Assuming that the right-handed neutrinos are produced thermally after in- 
flation results in a constraint on Mi [21], Mi < 10 Tr. This potential problem of thermal 
leptogenesis can be overcome e.g. in anomaly or gaugino mediated supersymmetry, as has 
been realized in p?51 139j. where the gravitino either decays before the initiation of nucle- 
osynthesis, or becomes the stable LSP. Since the gauge couplings decrease above a critical 
temperature T*, depending on the SUSY breaking scale, this mechanism leads to a relic 
gravitino density which is compatible with the WMAP results and which becomes indepen- 
dent of the reheating temperature for Tr > T*. In any case a small reheating temperature 
and Majorana mass Mi are desirable. We thus study in the following the values of the R 
matrix in (J48|) for which Mi acquires its minimum value. 

6 Constraints on the R matrix 

As outlined in Section 2, the R matrix ()11|) can be parametrized by six real parameters 
Xi, yi, i — 1,2, 3. With the Xi and yi given, Mi is fixed by the leptogenesis condition (J13J). 
In order not to violate the gravitino bound on the reheating temperature in mSUGRA 
models, we require Mi < 10 11 GeV. In addition, fhi has been demanded to lie in the interval 



Am 2 2 < fhi < yAm^, so that the fit for the washout factor Kf is appropriate. 



Another important condition for obtaining meaningful results on physical observables is 
the reliability of perturbative theory in the Higgs sector. We therefore require that the 
3 More stringent bounds on the reheating temperature can be found in the recent analysis |37| . 
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largest Yukawa coupling eigenvalue (Y u ) 3 meets the condition |(Yj,) 3 | 2 /47r < 0.3. The largest 
Yukawa coupling | [Y u ) 3 | is very sensitive to the heaviest Majorana mass M 3 and the pa- 
rameters Hi in the R matrix. This behavior can be easily understood from (jlUj) and the 
parametrization (jlljl which imply a y/M 3 dependence and an exponential dependence on 
Hi. Thus the perturbativity bound constrains yi < 0(1). In the subsequent numerical 
treatment an upper bound Ui < 1 has been imposed and remaining points with values of 
{Y v )l /4vr > 0.3 have been rejected. 

In the right-handed neutrino sector the following assumptions have been made. The R 
matrix parameters xi, xi and 23 are varied in their full range [0,27r] and the imaginary 
parts Hi are varied in the interval [10 -3 , 1]. For larger values of yi, the Higgs sector becomes 
strongly interacting. On the other hand, all three parameters yi have to be non-zero to 
make leptogenesis possible. Finally the hierarchical spectrum of the right-handed neutrinos 
has been generated by scattering M 3 logarithmically between 10 M± and 10 4 Mi and varying 
M2 between Mi and 0.1 M3. The condition M2 > ?>M% required for (42) to hold is met for 
the predominant majority of scatter points. 

The constraints on the R matrix imposed by the requirement of a minimal Majorana mass 
Mi for a given baryon asymmetry can be read off from (|43|) . It is obvious that M\ obtains 
its minimum values for Qm(Ru) ~ dte(Ru), and thus a small dte(Ru), since the imaginary 
parts are determined by the parameters While the contribution of R\\ in the sum in 
(}4*%j) is suppressed by the small value of m\ in hierarchical neutrino spectra, 3fte(i?i 2 ) and 
^te[Rxz) become small for £2,3 — mr with n G IN independently from the value of x±, as 
is obvious from pip. Thus, the remaining angle x\ cannot be constrained in this way. 
This is physically understandable, as x\ is related to the mixing between the right-handed 
neutrinos z/r 2 and i/r 3 , which plays no major role in the decays of v^. 
The allowed ranges of x<i and £3 are shown in Fig. Q for M\ = 10 10,11,12 GeV (left to 
right). All the other R matrix and light neutrino parameters are scattered in their ranges 
as outlined above and as given in Table ^ respectively. Each dot represents a viable data 
point which respects all the above selection rules. Fig. ^ also shows contours of constant 
Mi = 10 10 ' 11 ' 12 GeV as illustration for which these parameters are fixed at specific values. 
This behavior is also nicely illustrated in Fig. which shows the decrease of the allowed 
range of X2 with decreasing M\. Above a value of M\ pa 1.6 • 10 11 GeV, £2 can no longer be 
constrained and the full parameter range is allowed. The other R matrix angle x 3 behaves 
analogously. 
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Figure 1: Regions in the plane (x 2 , £3) consistent with the generation of the baryon asym- 
metry via leptogenesis and the gravitino bound M\ = 10 10,11,12 GeV, from left to right, 
respectively. The procedure yielding the scatter plots is described in the text. Also shown 
are the contours of constant Mi (^2,^3) = 10 10 ' 11 ' 12 GeV with all other parameters fixed 
(best fit values of oscillation parameters, mi = 0, vanishing CP phases, x± = and yi = 0.2, 
1 = 1,2,3). 

7 Charged Lepton Flavor Violation 



In the previous sections the requirement of successful leptogenesis has been used to constrain 
the lightest Majorana mass Mi and the parameters yt, x 2 $ of the R matrix. In the SUSY 
version of the seesaw model further observables are provided by the branching ratios of 
the radiative decays lj — > ^7. As has been shown in [2] in the framework of degenerate 
Majorana neutrinos, these branching ratios are particularly sensitive to the Majorana mass 
scale. In the present framework of hierarchical Majorana neutrinos, a study of Br(lj —>■ ^7) 
allows to obtain information about the heaviest Majorana mass M3. 

In the mass insertion approximation, the branching ratio Br (lj — > ^7) is schematically 
given by E] 



Br (lj -> l a ) 



a 3 tan 2 (3 mf 



(8m 



j' 



(49) 



where m denotes the typical mass scale of the sleptons in the loop. For our numerical 
calculations, we use the complete one-loop results for Br (lj — > ^7) [31 ITS]. 
Some general conclusions can be drawn from the dependence on the five relevant, unknown 
seesaw parameters (see (|15|h (|33|)) M 3 , 2r 12 and yi i2 , 



(8m 2 L )ij oc (YJLY U ) . . a M 3 | cos(x 2 + iy 2 ) 

\ / lj 



1 Mx f( j_ ■ ^ 



(50) 
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Figure 2: Region in the plane (X2, Mi) consistent with the generation of the baryon asym- 
metry via leptogenesis. The procedure yielding the scatter plot is described in the text. 



where 



fl2 
/23 



m 3 i \\ m 3 i 



(51) 
(52) 



These expressions can be contrasted with the present experimental limits |4TH I41j. 



-li 



Br(fji -> ei) < 1.2 • 10 
Br(r -> /i 7 ) < 6.8 ■ 10" 



as well as the expected sensitivity of the MEG experiment at PSI 
B-factories pFTj . 



Br(ji 
Br(r - 



pi) 



io- 13 , 

10~ 8 . 



(53) 
(54) 

and data from future 



(55) 
(56) 



From (|49I50|) it is immediately obvious that the branching ratios increase with the square of 
the heaviest Majorana mass M 3 . This behavior is illustrated in Fig.[3]for M% = 10 10,11 ' 12 GeV 
(left to right) and provides a possibility to determine M 3 from the combination of parameters 
M3 1 cos 62 1 2 — M 3 . The present bound on Br(fi — > erf) constrains M 3 already to be smaller 
than some 10 13 GeV. The sensitivity of the future MEG experiment at PSI would allow to 
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bound M 3 < O(10 12 ) GeV or determine the value of M 3 with an accuracy of a factor of 10. 
Due to the hierarchy M 3 /Mi > 10 used in our scattering, the minimum value for M 3 rises 
with Mi, pushing more and more points above the present limit on the branching ratio. 
The widening of the scatter range for Mi < 10 12 GeV (third plot) is largely due to the 
fact that %2 is no longer constrained around tot, as can be seen in the third plot of Fig. ^ 
Hence, Br(/i — > ej) oc | cos^l 2 can be strongly suppressed. 

The dependence Br(lj —> ^7) on the initial and final flavor is encoded in the functions /y, 
which show a sensitive variation with the value of x\. In particular, f\ 2 and f'23 given in 
(|51l52j) have minima at small yi and different values of xi, which imply a strong variation of 
the ratio Br(fi — > ej)/Br(r — > /ry) with x\. This ratio, which is rather independent of the 
mSUGRA scenario used ()33I49|) . is plotted in Fig. HJ The dependence on M 3 cos 2 #2 drops 
out and a measurement of both Br(fi — > ej) and Br(r — > ^7) would allow to determine 
Xi with an accuracy of about 0(1). It should be stressed, though, that present (|53I54|) and 
planned ()55I56|) experiments have a very restricted sensitivity, Br(fi — > ej)/Br(r —* yrf) < 
O(10 -3 ). If future B-factories would observe r — > /X7, the angle xi could also be restricted 
to values Xi ~ titc. 

8 Conclusions 

We have studied the minimal supersymmetric seesaw model by imposing constraints from 
light neutrino data, charged lepton flavor violation, successful leptogenesis and perturba- 
tivity in the neutrino sector. The parameters describing the light neutrino sector have been 
fixed to their experimental best fit values with an uncertainty anticipated from planned 
future experiments, and a very light neutrino spectrum, mj < m 2 < m 3 with mi < 0.03 eV 
has been assumed. 

The lightest right-handed neutrino mass Mi is fixed by the CP asymmetry e\ imposed to 
generate the observed value t]b = (6.3 ± 0.3) ■ 10~ 10 for the baryon asymmetry and the 
requirement of a small reheating temperature with Mi < 10 Tr < O(10 u ) GeV. The 
largest Majorana mass M 3 is constrained from above by the LFV process fi e^y. Thus 
the right-handed neutrino spectrum can be summarized as 

£>(10 n ) GeV ~ Mi < M 2 < M 3 < O(10 13 ) GeV, (57) 

in the given SUSY scenario SPSla. Alternative SUSY scenarios with sparticle masses in 
the range 100 GeV-1 TeV and low tan/3 are expected to yield similar results. 
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Figure 3: Branching ratio of jj — > cy as a function of the heaviest Majorana mass M3 
in SUSY scenario SPSla, for M\ = 10 10 ' 11 ' 12 GeV, respectively. The solid (dashed) line 
indicates the present (expected future) experimental sensitivity. For a description of the 
scatter procedure see text. 

In addition the successful scenario possesses an orthogonal R matrix encoding the mixing 
of the right-handed neutrinos, which is parametrized by the angles 9i = Xi + iyi with small 
but non-vanishing imaginary parts, being constrained as 



These bounds on y i: x 2 ,3 result from the requirement of perturbative Yukawa couplings and 
the condition of successful leptogenesis under the constraint of a small Mi. The remaining 
parameter x\ can be determined from the ratio Br(fi — > cy)/Br(r — > fij). 



When this paper has been accomplished the preprint jUj appeared, which comes to similar 




conclusions. One important difference between the two studies is that the authors of jlH] 



£2,3 — nrc, n G IN, 
< y t < 0(1). 



(58) 
(59) 



Note Added 
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7i/2 n 3/r/2 2n 

Figure 4: Ratio Br(/i — > e7)/Br(r — >■ /ij) as a function of x\ in SUSY scenario SPSla and 
for Mi = 10 11 GeV. For a description of the scatter procedure see text. Superimposed are 
two curves for which all parameters other than x± are fixed (best fit oscillation parameters, 
mi — eV, no CP phases and y 2 ,3 = 0). The solid (dashed) curve corresponds to y 1 = 0.01 
(2/i = 0.1). 

made further theoretical assumptions leading to more definite but less general results. 
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